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Hawking radiation is studied for arbitrary scalars, fermions and spin-1 bosons, using a tunneling 
approach, to every order in % but ignoring back-reaction effects, ft is shown that the additional 
quantum terms yield no new contribution to the Hawking temperature. Indeed, it is found that the 
limit of small fi in the standard quantum WKB approximation is replaced by the near-horizon limit 
in the gravitational WKB approach. 
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(3 ■ I. INTRODUCTION 

(N : 

\ Because it models Hawking radiation as the very physical process of quantum fields tunneling through a barrier, 

D ■ an approach called the tunneling method has recently gained popularity in the field of black hole thermodynamics. 

This method also provides multiple advantages vifhich go beyond merely providing intuition; indeed, because it 

, considers Hawking radiation to be a purely local phenomenon, it can be used to study spacetimes with multiple 

C\J horizons such as embedded black holes in deSitter spacetimes. Many spacetimes have thus been explored in this 

\ way: Kerr-Newman [l|, Black Rings Taub-NUT AdS black holes [5^, BTZjSJll, Vaidya dynamical 

' Q ' ' black holes |10| , Kerr-Godel [ll| , deSitter horizons [l^] , constant curvature black holes [iJl , as well as generic weakly 

qh| isolated horizons [l^ . 
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5—1 

{^JT^I Moreover, this method is especially powerful in that it allows quantum fields to be considered explicitly. As 
I— I such, the Hawking radiation of scalars to every order in Ti [l^ [l6|, spin-1/2 fermions to first order [l|, 0, IT7l - [l9| and 
later to every order [2(?|, higher-spin fermions [2l|,[22], and U{\) gauge bosons to second order [l^ has been inves- 
tigated. Expanding on these results, we will consider, to every order, the tunneling of scalars, fermions of any spin, 
and arbitrary gauge bosons from a generic near-horizon black hole metric. Some of the results we will present are not 
\^ , technically new; for example, the scalar field was calculated exactly in p^ . Nevertheless, we will include them not only 
T— I ■ for completeness, but also to provide a more thorough derivation and to interpret the results in a slightly different way. 

rn ; 

, The tunneling method comes in two fiavours. The first originates from the works of Volovik [2^ and later of 
T— I Kraus and Wilczek [23|, who analyzed the radiation process semi-classically by considering modes near the event 
, horizon; the idea was later generalized as a tunneling process by Parikh and Wilczek [25|. Near a Schwarzschild 
'■ horizon, radial null geodesies obey ^ = r = ±(1 — where the ± denotes the outgoing and incoming geodesies. 
J> ] The contribution to the imaginary part of the action comes from two part: a temporal contribution 2EImAt = Air ME 
from the discontinuity of the time coordinate at the horizon (26j . and a spacial contribution Im <fprdr, where we 
integrate along an infinitesimal complex path around the pole at the horizon. This closed path integral is to be 



. understood as a normalizing process which subtracts the infalling radiation from the outgoing one. Using the 
Hamilton equations of motion pr = we perform this integral to find Im/ = ImE f ^ — AttME. We then 

relate the tunneling rate to the action by F oc e~^™^ = e"®'^*^-^, thus finding the correct Hawking temperature of 

^ SttM- 

The second fiavour, which we will use throughout this paper, comes from the works of Padmanabhan and his 
collaborators [23]. This method was initially developed, and later formulated more algorithmically [7], as a 
means of studying the quantum tunneling of scalar particles through a gravitational barrier. It is rooted in the 
WKB approximation and consists of solving the Hamilton- Jacobi equations (thus earning it the nickname of 
Hamilton- Jacobi method) for a quantum field passing through an event horizon. For example, for a massless scalar 

field (f> = ae~ir^ , the equations of motion are — 0, of which the leading term is the Klein-Gordon equation 

^f,I^^'I + 0{h) = 0. This leads to / = ± / = ATiiME, yielding once more T = A goal of this paper is 

to show that this ?i — >■ approximation, which is also made for fermions and bosons, is unnecessary: since we are 
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only interested in near-horizon physics, we can use gu — ^ in place of ?i — to retrieve the correct Hawking tem- 
perature. This is somewhat unexpected since taking this limit is so common in semiclassical treatments of Hawking 
radiation; it means that we can consider this method to be a gravitational analog to the quantum WKB approximation. 

It is in that sense that our method calculates the temperature exactly to all orders. We will find that the 
action for a Klein-Gordon field obeys dri = '^JJT) ^^^ar the horizon. This implies two important points: first, since 
dtl is conserved, we can simply integrate this quantity to find the action without ever having assumed fi to be small: 
as such, our method is exact to every order in h. Second, this formula implies that if we expand the action / in powers 
of ?i as / = ^ h^Ii, then each li obeys this same equation: drU = ^777^ — "F / where the last equality defined 



the conserved quantity Ei = —dtii- This then means that drI ~ (l + J^i^i ^^^) ~ ^'-^0 (^1 + Si^i 



/(r) "^J /(r ^ 

E„ 



Besides h, there exists another parameter which is important to the problem of Hawking radiation: the ratio 
between the energy of the emitted particle and the mass of the black hole; this controls the amount of back- 
reaction on the black hole. Early in its development, the tunneling method was used to show that this back-reaction 
modified the thermal nature of the emitted radiation [2^. More recently, it was used to study correlations between 
emitted particles, and may provide a solution to the information puzzle to the lowest order [28l - [30l |. This area, 
however, is beyond the scope of the present paper. Thus, even though we will expand to every order in h, we will 
ignore back-reaction entirely. 

Our calculations will be done using a generic near-horizon line element in Schwarzschild-like coordinates: 

^ ~f{r)de + -^dr^ +dxl, (1) 

where /(r) vanishes at the horizon. This form is quite general and does not restrict us to spherically symmetric 
spacetimes. For example, in the near-horizon limit at fixed 9 ^ 6a, the Kerr metric can be written 

ds^ ^ -A(r, e)di^ + , + C(r, 9) \d<j) ~ D(r, 9)dtf + F{r, 9)d9'^ 

= -A,(ro, 9a){r - ro)dt^ + „ , , r + C(ro, 9o) [dcj) - ndtf , 

where = D{ro,9o) is the angular velocity of the black hole. This line element is of the form ([ij up to a 
redefinition of the r and coordinates. An unfortunate side-effect of using such a generic metric is a slight 
misdefinition of the energy, which we will illustrate here for the Kerr spacetime. From symmetry arguments, 
we know that the action is of the form / = —Et + J(j> + W{r,9). However, the redefinition of the (f> coordinate 
(j) = X + near the horizon means that the action is actually / = ~{E — ^lJ)t + Jx + W{r,9). Therefore, our 
energy actually corresponds to (E — J) in standard coordinates. Moreover, the temperature that we calculate is 
not redshifted: for an asymptotically flat spacetime, it represents the temperature measured by an observer at infinity. 

The purpose of this paper is therefore twofold. First, we will provide a generic treatment of bosons, which 
have so far only been studied in the Abelian case. Second, we will show that, ignoring back-reaction, terms of higher 
order in h do not modify the radiation process; in particular, we find that taking the near-horizon limit has the same 
effect as sending Ti 0. In the first three sections, we will analyze, respectively, Klein-Gordon, Rarita-Schwinger 
and Non- Abelian Yang-Mills fields radiating from the near-horizon metric ([T]) . Then, in Section |Vl we calculate the 
temperature associated with these particles, while, finally, in Section IVIl we discuss our conclusions and link our 
results with similar recent work in the field. 



II. SCALARS 

We begin by considering a massive scalar field cj) which we write (f> — . Although this form is based on the 
WKB approximation, we will not take the ?i — > limit which usually accompanies this approximation. The scalar 
field (f) obeys the Klein-Gordon equation: 

777^ 

= 9^"" (-9^9.0 + TIM) + -^^ 

f (3) 



3 



Upon reaching this stage, one commonly truncates the equation to leading order in h and takes the small-mass limit 
to retrieve the Hamilton- Jacobi equation d^Id^I = 0. This truncation, it turns out, is unnecessary, and we can 
analyze the situation to every order in Ti. Indeed, we begin by solving the above equation by assuming that 9^/ is 
real; the imaginary contribution to the action will come from integrating this divergent real quantity around a pole 
at the event horizon. We will ultimately find a solution for / which does solve the entire (complex) equation ([3]), thus 
justifying this assumption. Hence, taking the real part of equation fS]) and solving for dri , wc find 

a,,/ = ±Jz9l(^0j)2 + Zll^delf + ^{d^lY + (4) 

V 9^^ 9^^ V / grr ^ ^ ' gvr 

Because of the near-horizon symmetry, all of dfl, del and d^I are conserved quantities, and therefore finite. Moreover, 
g** will diverge, such that only the first term inside the square root will contribute. Thus, we find 



drl = ±^^dj. (5) 

As we mentioned earlier, this solution also solves the complex part of equation ([3]). Indeed, a simple calculation can 
show that, for drI defined by ([5]), we have 

g^'^d^dj = g'-'-a^/ = drI (ff"T;, + ff"rr,) = g^''T;,d.I, (6) 

since most Christoffel symbols vanish near the horizon. 

This result implies that / will obey the same equation even when we are not taking the limit of h going to 
zero. Instead, we take the limit of our metric approaching the event horizon; this is physically justified since that is 
where the tunneling takes place. Although more involved mathematically, we will find similar results for fermions 
and bosons in the following sections. 



III. FERMIONS 



It is intuitive that fermions must be emitted at the same temperature as scalar particles. Indeed, a fermion field 
of spin (n -I- i) is a tensor-valued spinor ^'^^..-^i^a which obeys the Rarita-Schwinger equations. Although these are 
commonly written, for a spin-3/2 field, as 

{e--P^^J^J,^p - ima'^^^) ^ ^a, (7) 
it is, for our purposes, much more enlightening to write them as 

H7^i^^-|-m)^'^,...^„, = 

where 

^M-5^-^V^[7.,7/3]; (9) 

this is the form in which they were originally studied [3lj |. In flat space, we can multiply by (i^'^D^ + m) to notice that 
each component of ^^i...^„a must obey the Klein-Gordon equation. Additional constraints then relate the components 
of with one another: the Dirac equation relates the Dirac indices a, while the other Rarita-Schwinger 

equations relate the higher-spin Lorentz indices fii ■ ■ ■ fin- Since we write the field as vE'^i...^„a — a^ii - ^i„ae^^ and 
are interested in the action, these extra relations play no role in the calculation of the Hawking temperature, and 
therefore fermions must be emitted at the same temperature as scalar particles. 

The Hawking temperature of fermions has already been calculated to leading order in ?i [l|, [13, 0, [Uj 
and, more recently, to every order [2^ for the massless case. We will perform the calculation for the massive 
arbitrary-spin case to every order in h, and will retrieve the results from j22| . albeit with an additional term which 
will not contribute to the Hawking temperature. The appearance of this term is due to our calculations being more 
thorough than previous ones, as we attempt to fill the gaps left behind by previous works. Moreover, our work has a 
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slightly different interpretation than that of [22f, as we will discuss in Section FVI] 
The Dirac equation implies that 



= -i-i^ d^a^,. 



(10) 



while the other Rarita-Schwinger equations will simply relate the various jii indices and will have no effect on the 
action; more details can be found in [2lj. We define the vierbein so that ejj^e'lrj^'^ = g^^; for the metric ([T]), this 

means e[j — \/\(f^5\. We also define the Dirac matrices 7^ = 6^7'', where the 7^ represent the flat-space 7 matrices 
in Majorana representation: 



1 

-1 



a' 
a' 



(11) 



where the tr' are the standard Pauli matrices. In particular, we have {7^,7''} — 2e^e;^77'"^ — 2g^'^ . Near the horizon, 
the metric only depends on the radial coordinate, which means that 9^ for ^ ^ r represents a Killing vector. This 
then implies that a^j...p^a can only be a function of r: the dependence of the fermion field on the other coordinates 
is restricted to a phase (such as e*^*). Therefore, '^^dfj,a^^...^^a = Ydra^-^...^^a — \fWl^ ^rO-i^i - i^^a. = 0, since 5'''' 
vanishes while a^-^...^^a remains finite. Hence, near the horizon, the Dirac equations become a system of equations 
linear in a,,,---u.,a- 



(12) 



= Det —rd^I + m+ -rg''^T'^^Aia,70] 



Reading this as a matrix equation in the Dirac indices, it is obvious that {-^l^d^I [70,7,9]) being 

invcrtible would imply aii^---p,^a = 0. Thus, we demand 



where we defined 



(13) 



—hm 





A 


B 





—hm 


C 


D 


-D 


B 


—hm 





C 


-A 





—hm 



Det 



{AD - BCf + 2m^{AD - BC) + m^ 



ih 



A = V^dti + Vg^dri + -r^/?^{9"gn'^^gu,r 



B 
C 

D ■ 



gOOdel + iy 
'^dgl-iy 



(14) 



ih 



As we approach the horizon, A and D diverge, such that the last terms do not contribute. We ultimately find AD = 0, 
which implies 




ih 



±E-:-{g''gn^/'gu,r 



(15) 



As discussed in [l7| . studying fermions provides us with insight which is absent from the scalar case: a direct 
meaning for the ± in (fT5|) . Indeed, consider the massless spin-1/2 case, where the fermion field is '^a = aa^^^ ■ The 
spin-up case corresponds to 



a 


/3 




(16) 
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where ^_|_ is the positive-spin eigenvector of ar- Then, combining equations ^TU\\ and ([T5|) . we find that either A = 
or B ~ 0. If A = 0, tta will be an eigenvector of 7^ with positive eigenvalue and therefore right-handed, whereas if 
B — 0, Ga will be left-handed. Thus, since they have the same spin, the two solutions of (fT5|) correspond to particles 
of opposite momenta: one is falling into the black hole whereas the other is outgoing. 



IV. BOSONS 



Although very little attention has been given to bosons using the tunneling rnethod, the emission of a C/(l) field 
from a generic black hole has recently been considered up to second order in h [22] . We will here expand on this result 
to find an exact formula for drl{r) for an arbitrary non-Abelian Yang-Mills theory. We therefore consider a vector 
field — a°e"B~^ obeying the equations of motion 



(17) 



= VF^,, 

where we defined F^^^ ^ d^A"^ - d^A^ + gf^^^A^^Al. Expanding this according to ^ oc aen^, wc find the expression 



= .9^ 



11 h n n 

cabc AC 



-gr^-Al[d^a'^--a''^d^I 



+ gr^Al d^al - -aid J 



- -aldxl - d.al + -aid J + gP'^A^Al 



( d^al - '-ald^I - d^al + -ald^I + g^'A^A^ 



'ihc Ab AC 



gr'^At I d^al - -ald^I - 9,a^ -f -a^9,/ + .g/^'^'^AM^ 



cde Ad AC 



We first simplify this expression by fixing the gauge: 



= V^A'^^ 



9a a° 



-a^d^I-Tt^al 



(18) 



(19) 



then, (ITSl) becomes 



O^g" 



dad^al - '-ay^d^I - '-d^ay^I 

i 1 
- dad^.a"^ + -daa^d^I ^ -^a^^dald^I 



-d a°-d I 



gr^Al d^a"^ - -a'^d^I 



- d^al - -ald^I - d^al + -aid J + g^'AlAl 



( dxal - -ald^I - d^al + gr'^A\Al 



gr'^At [ a„a^ - -a^9a/ ~ d^a^ + -a^d^I + g/^'^'^A^^A^ 



cde Ad AC 



(20) 
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We now focus on the ii = t equation. Setting the time derivatives of to zero and simphfying slightly, we find 
= 5"'' 



2h- 



1 Z 1 

-a'tdrdtl - drdrttf + -drafdrl H — Traldrldrl 
h n h 

abc Ah ' 



-afdrdrl + —drafdrl 



(21) 



+ gr'^Al i^drat - '-a^drl + '-a^dj + gr^'A^A'^ 
+ (•••), 

where the (■ ■ • ) refers to the 9 and sectors of the equation. We omitted those terms since they will not contribute 
to dri near the horizon, because they will remain finite while other terms, such as g'^*{dtl)^, diverge. Using the fact 
that the a° are normalized to be finite everywhere near the horizon (such that, for example, g^^a^^ — 0), and looking 
only at the real part of this equation (4^ , we find 



= g" 



+ ( 



(22) 



Since r*^ — ^g gu,r, then g 3> g^^^rt such that the two middle terms in the square brackets vanish. Our expression 
is therefore a quadratic polynomial in drI, which we will denote — ax^ + bx + c where a, b and ^ all go to zero at 

the same speed. Then, since - is finite while - diverges, we find x — ±\/ — . Hence: 



drI — ±\ 



-dj. 



(23) 



V. TEMPERATURE 



Calculating the Hawking temperature from drI has long been a contested issue, as many questions surrounding 
the covariance of the method have been raised. In particular, it appeared to yield different Hawking temperatures, 
differing by a factor of two [s?], depending on the coordinate system used. Now that the tunneling method has 
matured and become better understood, it is generally felt that we have a good handle on this issue. 



While many potential techniques have been proposed [34r|40[, we will here calculate the temperature using 
the approach recently summarized in [33| which, in particular, assumes that the tunneling rate follows a thermal 
distribution. The temperature gets two contributions: one from the integration over the radial coordinate and one 
from the discontinuity in the time coordinate: 

E 

" " lm{jdrl+- JdrI-+2EAt)' ^^^^ 

We begin by calculating the contribution from drI- In all cases, defining the energy as £' = —dfl, we have drI — 
jjp^ {±E + C) for some finite function C. It is clear, then, that 

E 
f [rH) 
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Next, we find the contribution from the discontinuity in the time coordinate. At, across the horizon. The metric ([T]) 
corresponds to an accelerated observer in flat space who follows the path 



tout = — smh(at) tin = — cosh(arj 

^ ^ (26) 

a^out = cosh(arj a^in = smh(arj, 

a a 

where the "in" and "out" subscripts refer to whether we are considering r < rn or r > rn , and where a = ^-j-^. 
Thus, as the horizon is crossed, we need t ^ t ^ so Ai = TJrj- Hence, 2EAt = f^Ty 8i.nd, from p4|. we get the 
Hawking temperature 

Th^^ (27) 

for every type of particle. This agrees with the temperature commonly found in the literature, which is usually 
calculated only to leading order in h. 

VI. DISCUSSION 

We've completed the study of spin-1 bosons, initiated in [2^ . by extending it to the non-Abelian case, by giving 
proper physical motivation for a number of terms dropping out, and by calculating all higher-order terms. Combined 
with previous results for scalars and fermions, this finally confirms that Hawking radiation is independent of the 
type of particle involved. Although our results show that bosons are emitted at the Hawking temperature regardless 
of the symmetries of the underlying theory, we needed to fix the gauge (V^A'' ~ 0) in order to perform the calculations. 

We've also shown that the ?i — >■ limit is unnecessary when using the tunneling method. Indeed, since the 
method is highly local in considering the emission of a field from a pole at the horizon, we are forced to take the 
limit r Th- It is therefore natural to take gu instead oi h 0, such that the tunneling method can truly 
be understood as the gravitational analog to the quantum WKB method. There are, however, some important 
drawbacks to this approach. First, by assuming that there is no back-reaction, we are drastically restricting quantum 
processes which may affect the radiation process. Second, one might expect to see grey-body corrections; the fact 



that these are missed by the tunneling approach makes this method suspect [4J] . 



It is also important to distinguish our results from those of Majhi et al. [Ta, 120, [241, '^^^ have recently cal- 
culated non-zero contributions from the Hawking temperature coming from higher-order terms in the tunneling 
method. This mismatch is simply a consequence of differing definitions of energy, which we can illustrate using 
the free scalar field. For this case, our calculations yield equation ([5|), which is also found in [l^. We define the 
energy as -E = —dtl, and therefore find no additional contributions to the Hawking temperature. On the other hand, 
^5] defines the energy as E = —dtio where Jq is the leading term in the action / = '^fi'Ii, this then produces 
higher-order corrections. We are not the first to question these corrections; indeed, [il!,!!^] have pointed out that they 
are caused by an odd definition of the field's energy and concluded that the Hawking temperature is not modified by 
higher-order terms. Moreover, in previous works on this topic, terms with no explicit dependence on the action (such 
as the last term in our equation (jlSp ) are automatically set to zero; we've filled this gap by providing the necessary 
justification as to how each term cannot contribute to the Hawking temperature. 

In conclusion, we have calculated the temperature associated with the emission rate of every type of scalars, 
fermions, and spin-1 bosons from a generic black hole spacetime and to every order in T% in the tunneling approach, 
as long as there is no back-reaction. 
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